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ABSTRACT 

We use the GALPROP code for cosmic -ray (CR) propagation to calculate the broad-band luminosity spec- 
trum of the Milky Way related to CR propagation and interactions in the interstellar medium. This includes 
7-ray emission from the production and subsequent decay of neutral pions (7r°), bremsstrahlung, and inverse 
Compton scattering, and synchrotron radiation. The Galaxy is found to be nearly a CR electron calorimeter, 
but only if 7-ray emitting processes are taken into account. Synchrotron radiation alone accounts for only one 
third of the total electron energy losses with ~ 10 — 20% of the total synchrotron emission from secondary 
CR electrons and positrons. The relationship between far-infrared and radio luminosity that we find from our 
models is consistent with that found for galaxies in general. The results will be useful for understanding the 
connection between diffuse emissions from radio through 7-rays in "normal" (non-AGN dominated) galaxies, 
as well as for estimating the broad-band extragalactic diffuse background from these kinds of galaxies. 

Subject headings: Galaxy: general — gamma rays: galaxies — gamma rays: general — radio continuum: 
galaxies — radiation mechanisms: non-thermal — ISM: cosmic rays 



1. INTRODUCTION 

Cosmic rays (CRs) fill up the entire volume of galax- 
ies, providing an important source of heating and ionisation 
of the interstellar medium (ISM), and may play a signifi- 
cant role in the regulation of star formation and evolution 
of galaxies (lFerrierell200l l ICoxll200l ISocrates et aD l2008l 
ISironi & Socratesl 120101) . Diffuse emissions from radio to 
high-energy 7-rays (>100 MeV) arising from various inter- 
actions between CRs and the ISM, interstellar radiation field 
(ISRF), and magnetic field, are currently the best way to trace 
the intensities and spectra of CRs in the Milky Way (MW) 
and other galaxies. Gamma rays are particularly useful in this 
respect since this energy range gives access to the dominant 
hadronic component in CRs via the observation of 7r°-decay 
radiation produced by CR nuclei inelastically colliding with 
the interstellar gas. Understanding the global energy bud- 
get of processes related to the injection and propagation of 
CRs, and how the energy is distributed across the electromag- 
netic sp ectrum, is essential t o interpret the radio/far-infrared 
relatio n (Ide Jong et al.lll985b iHelou et all 19851: iMurphy etai] 
2006), galactic calorimetry (e.g JVolk 11989b. and predictions 
of extragalactic ba ckgrounds (e.g., iThompson et al.1 120071: 
IMurphy et ai1l2008b . and for many other studies. 

The MW is the best studied non-AGN dominated star- 
forming galaxy, and the only galaxy that direct measurements 
of CR intensities and spectra are available. However, because 
of our position inside, the derivation of global properties is 
not straightforward and requires detailed models of the spa- 
tial distribution of the emission. Nevertheless, constructing a 
model for the global properties of such a galaxy is tractable 
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with the variety of data available. 

In this paper, we calculate the injected CR power and cor- 
responding broad-band luminosity spectrum from radio to 7- 
rays for CR propagation models consistent with current CR, 
radio, and 7-ray data. Earlier estimates focussing on the MW 
7-ray luminosity only (e.g.. IStrong et al.l 1200 0) were based 
on modelling EGRET data. The launch of the Large Area 
Telescope (LAT) on the Fermi Gamma-ray Space Telescope 
(hereafter, Fermi-LAI), has provided a wealth of new 7-ray 
data up to, and beyond, 100 GeV. Analysis of the ferm/-LAT 
data has not confirm ed the anomalous "EGRET GeV-excess" 
dAbdo et al. 2009a) and has led to an improved model of CR 
propagation and diffuse 7-ray emission thus enabling a better 
estimate of the CR and 7-ray luminosities of the Galaxy. 

2. THE MODELS 

We use a model based on the GALPROP code0 for the 
CR-related processes in the Galaxy that has been adjusted 
to fit many types of data including direct measurements of 
CRs, non-thermal radio emission, hard X-rays, and 7-rays. 
Studies of the diffuse Galactic emission and CRs using this 
code prior to the laun c h of t h e Fermi-LAT can be found i n 
I Moskalenko & Strong! d 1998b. IStrong & Moskalenkol d 1998b 
Istrong et al.l (12000b. IMoskalenko et all (12002b. IStrong et all 
(2004a), Strong et al.M2004bl) . and iPorter et al.1 d2008b . An 
extensive revie w of CR propagation , models, data, and litera- 
ture is given by IStrong et al.l ( 12007b . 

We calculate diffuse emissions from radio to 7-rays, pro- 
duced by CR protons, helium, and electrons/positrons. The 
contributions by discrete sources and line emissions (e.g., 
511-keV annihilation radiation) are not included. Since the 
propagation parameters are not uniquely determined from the 
observations we consider an illustrative range of parameters 
for diffusive-reacceleration (DR) and plain diffusion (PD) 
propagation models. The main uncertainty we consider is the 
CR confinement volume and associated propagation model 
parameters, e.g., the diffusion coefficient. Other effects to 
consider would include different propagation modes, such as 

6 http://galprop.stanford.edu 
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convection, the distribution of CR sources, etc. However, 
these are outside the scope of the current paper and we de- 
fer these to future work. 

The output of a typical GALPROP run includes CR dis- 
tributions, 7-ray/synchrotron skymaps and emissivity distri- 
butions. The luminosities of the different components are 
computed by integrating the respective emissivities over the 
total Galactic volume: for IC and synchrotron radiation the 
volume emissivity is directly calculated, while for ir° -decay 
and bremsstrahlung the emissivity is calculated per hydrogen 
atom so we weight the calculated volume emissivities by the 
distributions of atomic, molecu lar, and ionised hydrogen in- 
cluded in the GALPROP code dStrong & Moskalenkoll 19981; 
iMoskalenko et alj|2002trstrong et al.ll2004bl) . 

We obtain the input CR luminosity for each of our models 
by integrating the model CR source spectra and spatial dis- 
tribution over the Galaxy. The injected luminosity, instead 
of emergent, is the relevant quantity because we are inter- 
ested in how the CR power is transferred to the electromag- 
netic channels. In addition to CR protons and helium, and 
primary electrons, we also calculate the luminosities of sec- 
ondary positrons and electrons resulting from CR-gas interac- 
tions in the ISM. The diffuse emissions from these secondary 
CR species are also calculated and included in the total syn- 
chrotron, bremsstrahlung, and IC spectra. 

We use DR and PD propagation models where the spa- 
tial diffusion coefficient D xx and its momentum dependence, 
together with the diffusive reaccelation characterised by an 
Alfven speed, va, assuming a Kolmogorov spectrum of inter- 
stellar turbulence (if used), and the size of the CR confinement 
volume, Zh, are obtained by fitting the CR second ary/primary 
ratios for B/C and 10 Be/ 9 Be; for details see Stro ng et all 
(2007). The CR source injection spectra are taken as broken 
power-laws in momentum, with different parameters for nu- 
clei and primary electrons. These are chosen to reproduce the 
directly observed CR spectra after propagation for the adopted 
models. 

The calculations are based on model pa rameters that repro- 
duce the Fermi ~LAY electron spec tr um dAbdo et alj [2009a) 
and 7-ray data (Ab do et alJ l2009a b. 2 010allbl) . The spatial 
distribution of the CR sources is based on pulsars as tracer s 
of supernova remnants (SNRs) used in Stro ng et alJ (|2004b), 
but constant beyond the solar circle as indicated by analysis of 
the Fer w/-LAT data for the 2 nd Galactic quadrant (I Abdo et afl 
l2010ah . 

The synchrotron calculation uses the same CR electron 
model and a magnetic field model that reproduces radio sur- 
veys at frequencies from 100 MHz to 23 GHz (Strong et al. 
in preparation). Diffusive-reacceleration and PD propagation 
models are used for three halo sizes, Zh = 2, 4, and 10 kpc, 
respectively, with corresponding self-consistently derived dif- 
fusion coefficients D xx {p) — Do(3{p/ po) S , where (3 = v/c 
is the dimensionless particle velocity, p is the particle rigid- 
ity, and Do, po, S are constants. These halo sizes cover the 
range consistent with available CR data for B/C and 10 Be/ 9 Be 
(Enge lmanneTall [19901 lYanasak et all HOQlt Ide Nolfo et ail 
2006, and references therein). The parameters and values for 
the models are summarised in Table Q] 

The optical to far infrared (FIR) luminosity spectrum is de- 
rived from the model of the ISRF used for the propagation 
calculations (the " maximum metallicity gradient" model from 
iPorter et al. 2008). The emergent luminosity for the ISRF is 
computed by surface integration over a region large enough 



to encompass the total flux from the stellar luminosity distri- 
bution and the starlight reprocessed by dust (^30 kpc radius 
about the Galactic centre). In the present work, the input bolo- 
metric stellar luminosity is ^4 x 10 10 L© apportioned across 
the stellar components boxy bulge/thin disc/thick disc/halo 
with fractions ~0. 1/0.7/0. 1/0.1, and ^20% reprocessed by 
dust and emitted in the infrared. The uncertainties related to 
the distribution of the ISRF interior to the integration bound- 
ary include the relative luminosities of the bulge component 
and the disk, the metallicity gradient, and other details. How- 
ever, these uncertainties are dominated by the overall uncer- 
tainty in the input stellar luminosity. A higher input stellar 
luminosity will increase the CR electron/positron losses via 
IC scattering and hence the overall output in 7-rays, requiring 
a higher injected CR power and increasing the overall calori- 
metric efficiency (see below). Estimates available in the lit- 
erature illustrating the r ange for the MW stellar luminosity 
are, e.g., 6.7 x 10 10 L e dKent etalJll99lh and 2.3 x 10 10 L Q 
dFreudenreichll998l) . 

3. RESULTS AND DISCUSSION 

Figure Q] shows the broad-band luminosity spectrum of the 
Galaxy, including the input luminosity for CRs for a 4 kpc 
halo for a DR and PD model, respectively^. Fig.|2]illustrates 
the detailed energy budget for the DR model with 4 kpc halo 
size. Figure [3] shows the broad-band luminosities for DR and 
PD models for the 3 halo sizes, and Table [2] summarises the 
spectrally-integrated luminosity for the various components 
for each of our models. The peak injected CR luminosities 
differ by a factor ~ 2 for different propagation modes, but the 
total injected CR luminosities are close, only differing at the 
~few precent level. The injected CR spectra differ between 
DR and PD models to compensate the momentum dependence 
of the diffusion coefficient in each propagation model so that 
the local CR spectra are reproduced. Independent of propaga- 
tion mode, the relative decrease in the injected CR proton and 
helium luminosities is ~ 10% for halo sizes 2 to 10 kpc. For 
smaller halo sizes, the CRs escape quicker requiring more in- 
jected power to maintain the local CR spectrum. In addition, 
for larger halo sizes CR sources located at further distances 
can contribute to the local spectrum, which is our normalisa- 
tion condition, hence less power is required. In contrast to 
nuclei, the injected primary CR electron luminosity increases 
with Zh, reflecting the increased input power of these particles 
required to counter the energy losses from the larger confine- 
ment region and escape time. 

The CR nuclei luminosities can be directly compared with 
the well-known appr oximate estimate fr om CR "grammage" 
as described by, e.g.. lDogiel et al.1 (120021) . which shows these 
can be estimated from the observed secondary/primary ratios. 
The values in Tableware within a factor 2 of the ^10 41 erg 
s -1 derived in that paper, providing a consistency check on 
our more detailed modelling. Taking the CR lumi nosities and 
adopting a supernova rate of (1.9 ± l.l)/century (Die hl et alJ 
120061) we obtain a CR energy input per SNR of (0.3- 1) x 10 50 
erg, comparable to standard estimates of a 5-10% efficiency 
for CR acceleration in SNRs. 

For the DR model, the injected CR lepton luminosity is 
close to the total of synchrotron, IC, and bremsstrahlung lu- 
minosities for the 10 kpc halo (Model 3), showing that the 
Galaxy is approximately a lepton calorimeter (^79% effi- 

7 Spectra for all models and components will be made available in numer- 
ical form via the German Virtual Observatory, http://www.g-vo.org 
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ciency) in this case, with the bulk of the energy lost via IC 
emission (only about one third is lost via synchrotron emis- 
sion). Smaller halos still constitute fair calorimeters but are 
less efficient: ~51% for the case of a 2 kpc halo with ~17% 
of the energy loss due to synchrotron radiation, and ^59% for 
the case of a 4 kpc halo with ~19% of the energy loss due to 
synchrotron radiation. For PD models the calorimetric effi- 
ciencies are decreased, but the same trend of higher efficiency 
with increasing halo size is evident. 

Our calculations show that calorimetry holds for CR lep- 
tons in MW-like galaxies independent of propagation mode 
provided the size of the CR confinement volume is large 
enough to allow for sufficient cooling and if 7-ray produc- 
tion is included. A large fraction of the IC luminosity is at 
energies below 100 MeV, extending to hard X-rays, and in- 
cludes a significant secondary electron /positron component 
(IStrong et al.l l2004at iPorter et al.1 120081) . The contribution 
to the synchrotron emission from secondary electrons and 
positrons is ~ 10 — 20% for our Galaxy-like models, being 
largest for small halo sizes. For DR models, the secondary 
lepton emission is 20%, 16%, and 16% for 2, 4, and 10 kpc, 
respectively. For PD models, the secondary lepton emission 
is 12%, 10%, and 10% for 2, 4, and 10 kpc respectively. 

Using the standard definition of the FIR/radio relation, q = 
log 10 (S , Fm/3.75 x 10 12 Hz) - log 1Q (Si. 4 GH*)> where Spin, 
(W m" 2 ) is the FIR flux from 42.5-122.5 /zm a nd 5i. 4 qhz 
(W m ~ 2 Hz -1 ) is the radio flux at 1.4 GHz dHelou et alJ 
[19851) . we calculate this ratio for our luminosity models (see 
Table [2] under "radio-FIR relation"). The range of values 
2.26 — 2.69 is consistent with that for the correlation over 
normal and starburst galaxies: 2.34 ± 0.26 (lYun et al. 11200 ll) 
and near to the value for M33 (2.50) given in that paper. Note, 
here we have used the FIR luminosity from the ISRF model 
described earlier, 1.9 x 10 43 erg s _1 (see Table 0. How- 
ever, the FIR luminosity of the Galaxy is uncertain by at least 
a factor 2 and t his w ill affect q. Using COBE/DIRBE data, 
ISodroski etaD (119971) give a total IR luminosity 4.1 x 10 43 
erg s _1 and IR luminosity > 40 /im, IR>4o M m> of 2.8 x 10 43 
erg s _1 . This yields a FIR luminosity of 1.4 x 10 43 erg s _1 
from FIR/I R^Qmn = 0.5 1 for th e ISRF used in this paper. 
Meanwhile, Paladin i~et al.l (120071) give a much higher value 
for IR>60pm, 1.1 x 10 44 erg _1 , yielding a FIR luminosity 
6.6 x 10 43 erg s _1 (Paladini, private communication). This 
higher value for the FIR luminosity would increase all q val- 
ues in Table|2]by 0.53. 

The luminosity spectra for our models contains significant 
IC 7-ray emission, the contribution at high energies depend- 
ing on details of the propagation model. Above ~10 GeV the 
contribution by IC is non-negligible compared to 7r°-decay. 
Estimates of the isotropic background using models for the 
diffuse 7-ray emission of normal star-forming galaxies invok- 



ing only 7r - decay as the dominant m ode for high-energy pro- 
duction, e.g., Ando & Pavlidou ( 2009), and ignoring details of 
the propagation will underestimate the contribution per galaxy 
by these classes of objects. 

The present estimates are intended to be illustrative, using a 
particular set of models; the level of uncertainty is suggested 
by the range covered by the models which are far from ex- 
haustive. In any case the general conclusions from this work 
are robust. 
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for useful comments. T. A. P. acknowledges support from 
NASA Grant No. NNX10AE78G. I. V. M. acknowledges sup- 
port from NASA Grant No. NNX09 AC 1 5G. 
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Table 1 

GALPROP model parameters 







Diffusive Reacceleration 






Plain Diffusion 






Model 1 


Model 2 


Model 3 


Model 1 


Model 2 


Model 3 


Model parameter/GALDEF ID 


ZU2LMS 


z04LMS 


zlOLMS 


z02LMPDS 


z04LMPDS 


zlOLMPDS 


Halo height, kpc 


2 


4 


10 


2 


4 


10 


Galaxy radius, kpc 


20 


20 


20 


20 


20 


20 


Diffusion coefficient" Do 


2.9 


5.8 


10.0 


1.8 


3.4 


6.0 


Diffusion coefficient" S 


0.33 


n 11 


n ii 


D.J 


U.j 


0.5 


Reacceleration va, km s — 1 


30 


30 


30 








CR sources distribution 








Pulsars (Lorimer 2004) b 






Magnetic field strength 






B = 








Injection spectrum (nuclei) 














Index below break 


1.98 


1.98 


1.98 


1.80 


1.80 


1.80 


Index above break 


2.42 


2.42 


2.42 


2.25 


2.25 


2.25 


Break energy, GeV 


9 


9 


9 


9 


9 


9 


Normalisation energy, GeV 


100 


100 


100 


100 


100 


100 


Proton normalisation intensity d,c 


5 


5 


5 


5 


5 


5 


Injection spectrum (primary electrons) 












Index below break 


1.60 


1.60 


1.60 


1.80 


1.80 


1.80 


Index above break 


2.42 


2.42 


2.42 


2.25 


2.25 


2.25 


Break energy, GeV 


4 


4 


4 


4 


4 


4 


Normalisation energy, GeV 


34.5 


34.5 


34.5 


34.5 


34.5 


34.5 


Normalisation intensity* 1 


0.32 


0.32 


0.32 


0.32 


0.32 


0.32 



a D xx = 10 28 /3Do(p/ po)* 5 cm 2 s 1 , po = 4 GV, /3 = v/c, constant below po for plain diffusion model. 

b Modified to have the the value at R = 10 kpc across the range R = 10 — 15 kpc and to be zero beyond 15 kpc. 

c Ro = 8.5 kpc, Rb = 50 kpc, zb = 3 kpc. 

d 10 — 9 cm -2 s — 1 sr _1 MeV -1 (nucleon -1 for nuclei). 

c Normalisation energy 100 GeV for protons corresponds to a rigidity ~ 100 GV. Helium normalisation is 0.068 relative to protons at the same rigidity. 
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Table 2 










Luminosity of the Galaxy for various processes, 10 38 erg s 












Diffusive Reacceleration 




rlain Diffusion 






MnHpl 1 


Model 2 


Model 3 


Model 1 


Model 9 


Model 3 


Cosmic rays (0.1-100 GeV): 


805 


790 


698 


780 


723 


660 


Protons 


737 


724 


633 


718 


662 


601 


Helium 


56 


55 


48 


52.4 


48.3 


43.9 


Leptons 


12.2 


14.5 


16.9 


1 n da 


1 9 Q 1 

iz.y i 


1 J. 1 J 


Primary e 


8.8 


11.1 


13.4 


£ AS 
o.OJ 


it/, j 


IZ./ 


Secondary e~ 


0.78 


0.77 


0.83 


0.63 


0.04 


n as 
O.OJ 


Secondary e~*~ 


2.6 


2.6 


2.7 


J ~7A 
l./O 


J TV 


J 0/1 

l.oO 


7-rays (0.01-100 MeV): 


2.32 


3.34 


6.22 


1.47 


2.20 


3.50 


7r° -decay 


0.24 


0.23 


0.23 


0.14 


0.13 


0.13 


Inverse Compton 


1.80 


2.81 


5.63 


1.09 


1.79 


3.04 


Primary e~ 


1.31 


2.20 


4.41 


0.91 


7.54 


2.67 


Secondary e 


0.49 


0.61 


1.22 


0.18 


0.25 


n ?*7 
U.J / 


Bremsstrahlung 


0.27 


0.30 


0.36 


U.Z*t 


U.Zo 


U. jj 


Primary e~ 


nil 
0.11 


0.15 


0.19 


0.15 




n i/t 
U.Z4 


Secondary e 


0.16 


0.15 


0.17 


n no 
O.Oy 


/i no 
O.Oy 


n no 
O.Oy 


7-RAYS (0.1-100 GEV): 


8.86 


9.12 


10.3 


6.72 


7.18 


7.79 


7r°-decay 


6.75 


6.46 


6.59 


4.99 


4.90 


4.79 


Inverse Compton 


l.Zj 


1.76 


2.59 


1.27 


1.77 


2.43 


Primary e~ 


1.1 J 


1.66 


2.43 


1.20 


1.68 


2.33 


Secondary e^ 


If) 


0.10 


0.16 


u.u/ 


n no 
U.Uy 


n m 


DICIIlSt»LI dillUIlg 


0.87 




1 oc 




U.J 1 


V.J 1 


Primary e~ 


0.51 


U.JO 


L/. /*f 


0.33 


0.39 


0.45 


Secondary e- 11 


0.36 


0.30 




0.13 


0.22 


0.12 


KAU1U ^U.UU1-1 UU OrlZ. ) . 


2.07 


HA 
L. /D 


J. /Z 


1 Ait 


z.io 


Z.jU 


Primary e — 


1.65 


2.31 


3.10 


I.3U 


1 Q 1 
1.01 


z.zo 


kjCL/Uiiucu y c 


0.42 


U.4j 


U.Oz 


18 


0.22 


0.24 


Conversion efficiencies: 














f loylv^lx ICULUHA 


0.34 


0.40 


0.57 


28 


0.34 


42 




0.17 


0.19 


0.22 


0.13 


0.16 


0.17 


Lepton calorimetric efficiency c 


0.51 


0.59 


0.79 


0.41 


0.49 


0.59 


HpTTr A T PTR (0 ^ ^flOfl TH7V 






24.8 X 10 43 erg s" 


-l 






upticai 






19.6 x 10 43 erg s" 


-l 






local IK 






5.2 X 10 43 erg s" 


i 






FIR f 






1.9 X 10 43 ergs" 


i 






RADIO-FIR RELATION: 














Synchrotron (1.4 GHz)« 


1.68 


2.18 


2.85 


1.06 


1.49 


1.82 


9FIR 


2.49 


2.38 


2.26 


2.69 


2.54 


2.45 



a (IC + bremsstrahlung)/(total leptons). 
k (Synchrotron)/(total leptons). 

c (IC + bremsstrahlung + synchrotron)/(total leptons). 

d 0.1 - 8 (im; 38 - 3000 THz. 

e 8- 1000 ^tm;0.3 - 38 THz. 

f 42.5 - 122.5 Jim; 2.45 - 7.06 THz. 

g 10 28 erg S - 1 Hz -1 . 
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.galdef ID 54_z04LMPDS 
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log (Energy, MeV) 
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Figure 1. Global CR-induced luminosity spectra of the MW for DR propagation model (left panel) and PD propagation model (right panel) with z h = 4 kpc. 
Line styles: ISRF, including optical and infrared scaled by factor 10~ 4 (magenta solid) and components for model 2 - Cosmic rays (dotted lines), protons (red), 
helium (blue), primary electrons (green), secondary electrons (cyan), secondary positrons (magenta); CR-induced diffuse emissions (solid lines), IC (green), 
bremsstrahlung (cyan), 7r°-decay (red), synchrotron (black, left side of figure), total (black, right side of figure). 
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Figure 2. The luminosity budget of the MW for DR propagation model with = 4 kpc. The percentage figures are shown with respect to the total injected 
luminosity in CRs, 7.9 X 10 40 erg s — 1 . The percentages in brackets show the values relative to the luminosity of their respective lepton populations (primary 
electrons, secondary electrons/positrons). 



Luminosity and energy budget of the Milky Way 



7 




log (Energy, MeV) 




log (Energy, MeV) 



Figure 3. Global CR-induced luminosity spectra of the MW for range of halo sizes for DR propagation model (left panel) and PD propagation model (right 
panel). Line styles: Zh = 2 kpc (solid), z^ = 4 kpc (dotted), and = 10 kpc (dot-dashed). Components as in Fig.fT] 



